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Abstract: We have measured the charged kaon hfetime using a sample of 15 x 10^ 
tagged kaon decays. Charged kaons were produced in pairs at the DA$NE c/i-factory, 
e^e~ — > ^ K^K" . The decay of a was tagged by the production of a and 
viceversa. The hfetime was obtained, for both charges, from independent measurements of 
the decay time and decay length distributions. From fits to the four distributions we find 
r = (12.347±0.030) ns. 
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1. Introduction 

While the lifetime of the Ks is well measured and measurements of the lifetime of the 
Kl [1], [2] have been recently performed, the most precise measurement of the lifetime 
of the charged kaon dates back to 1971: t± = (12.380 it 0.016) ns [3]. At the time, the 
agreement with previous measurements was poor and more recent measurements did not 
improve the agreement. The 2006 PDG [4] gives t± = (12.385 it 0.025) ns corresponding 
to an accuracy of 0.2%. However the set of measurements upon which the above value is 
based is not self consistent. The probability that the spread of the values used by PDG be 
due to the statistical fluctuations is very low, 0.17%. Averaging an inconsistent set is not 
a valid procedure. The PDG in fact enlarges the errors by a factor of 2.1. It is important 
therefore to confirm the value of t±. We report in the following on new measurements of 
T±. The statistical error on the lifetime obtained fitting a decay curve over a time interval 
At is given by^: 

5t I - 1 , , 

— = —= X , (1.1) 
T y/N ^l + e2^-e^(2 + r2) 

^this formula is just the standard maximum likelihood estimator error evaluation for an exponential 
probability distribution function integrated over a finite time interval At 
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where T = At/r is the time interval in ii' -hfetime units. With T ~ 2 and number of 
events N a A x 10®, the best statistical accuracy reachable in our case is 0.1%, if no other 
sources of statistical error were present. We have developed two different methods, one 
employing the reconstruction of the kaon path length and the other measuring directly the 
kaon decay time. 

2. The KLOE detector 

Data were collected with the KLOE detector at DA$NE [5], the Prascati e+e" collider 
operated at a center of mass energy W = ~1020 MeV. Equal-energy positron and 
electron beams collide with a crossing angle of vr — 25 mrad, which results in a small 
transverse momentum ( ~ 12.5 MeV) of the (p mesons produced. The beam spot has 
cjj; ~ 2 mm Gy ~ 0.02 mm and o"^ ~ 30 mm. The KLOE detector is inserted in a 0.52 T 
magnetic field. It consists of a large cylindrical drift chamber (DC), surrounded by a fine 
sampling lead-scintillating fibers calorimeter (EMC). The DC [6], 4 m diameter and 3.3 m 
long, has full stereo geometry and operates with a gas mixture of 90% Helium and 10% 
Isobutanc. Momentum resolution for tracks with large PT is t{p±)/pt ~ 0.4%. Spatial 
resolution is arct, ~ 150/xm and cr^ ~ 2 mm. Vertices are reconstructed with a 3D accuracy 
of ~3 mm. The DC inner walls arc made in carbonium fibers with thickness of ~1 mm. 
The beam pipe walls are made of AlBeMet, an alloy of beryl-aluminum 60%-40%, with 
thickness of ~0.5 mm. The EMC [7], divided into a barrel and two endcaps, for a total 
of 88 modules, covers 98% of the solid angle. Arrival times of particles, 3D positions and 
the energy deposits are obtained from readout of signals on both ends of the module, with 
a granularity of 4.4 x 4.4 cm^, for a total of 2240 cells arranged in five layers. Cells 
close in time and space are grouped into a calorimeter cluster. Resolution on energy and 
time measurement are as/E = 5.7%,/ ^/E(GeV) and or = 57ps/v/E(GeV) ® 100 ps. The 
trigger [8], used for this analysis, is based on the coincidence of at least two local energy 
deposits in the EMC, above a threshold of 50 MeV in the barrel and 150 MeV in the 
end caps. Cosmic rays muons crossing the detector are vetoed. Since the trigger time 
formation is larger than the interbunch time ('^ 2.7 ns), the trigger operates in continuous 
mode and its signal is synchronized with the DA$NE accelerator Radio Frequency. The 
trigger Tq for an event is thus displaced in time with respect to the correct crossing by 
an integer multiple of the interbunch time which is different event by event. The correct 
bunch crossing is found after the event reconstruction. For this analysis a sample 3x10^ of 
— K'^K~, events generated with the KLOE MonteCarlo (MC) [9], has been used. 

3. The tag procedure 

0-mesons decay, in their rest frame, into anti-collinear KK pairs. In the laboratory this re- 
mains approximately true because of the small ^-meson p. The detection of a K~^{K~) tags 
the presence of a K^{K^) of given momentum and direction. The decay products of the 
K^K^ pair define two spatially separated regions called the tag and the signal hemispheres. 
Identified decays tag a sample of known total count N ■ This procedure is a unique 
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feature of a 0- factory and allows to measure absolute branching ratios. For this analysis 
charged kaons are tagged using the — >■ /U^^i^]^ {^1x2) decay. This decay corresponds 
to about 63% of the charged kaon decay width [4] and since BR{(j) — K'^K~) 2± 49% [4] 
and a{e'^e~ (f)) ^ 3 /xb, there are about 1.1 x lO^K'^K" events/pb~^. The iC^2 de- 
cay is clearly identified since it exhibits a peak, reconstucted with a resolution of about 
1 MeV, at ~ 236 MeV in the momentum spectrum of the secondary tracks in the kaon 
rest frame. In order to minimize biases to trigger efficiency, we ensure that the tagging 
kaon did provide the trigger of the event by requiring a cluster in the EMC (associated to 
the muon track) which satisfies the trigger conditions. Hereafter these events are called 
self-triggered tags. We find ATgeift^g tag ~ 1-2 x 10^ per pb~^. The MC simulation shows 
that the contamination of the selected cf) — ^ K^K~ sample is negligible. The K^^z tag 
allows a precise determination of the correct bunch crossing of the event using the muon 
and kaon track lengths and the momenta measured in the DC and the arrival time of the 
muon in the EMC. The tagging efficiency is almost independent on the proper decay time 
of the signal kaon. The small residual correlation has been evaluated with the Monte Carlo 
and checked with data using doubly tagged events, the events in which both the and 
the K" are reconstructed and tag the event. 

4. Signal selection 

The measurement is performed using data from an integrated luminosity C= 210 pb~^ 
collected at the (j) peak. The average c/)— meson momentum and the coordinates of the 
e^e^ interaction point (IP) are measured run by run with Bhabha scattering events. K^2 
tags of both charges have been used. We developed two analysis methods: the kaon decay 
length and the kaon decay time. The two methods have comparable precision and different 
systematics; this allows a useful cross-check of the results. We use a coordinate system 
where the x-axis points to the center of the collider, the z-axis bisects the two beam lines 
and the y-axis is vertical. For both methods the kaon decay vertex position (x, y, z) is 
searched for in a fiducial volume (FV) defined by: 



where R = yx^ + y^. 

4.1 Kaon decay length method 

The measurement of the charged kaon decay length requires the reconstruction of the kaon 
decay vertex using only DC-information. We use any kaon decay mode. The evalua- 
tion of the vertex reconstruction efficiency uses a data control sample obtained by EMC- 
information only. Given a charged kaon fulfilling the self-triggering tag requirements, our 
signal is given by the opposite charged kaon decaying in the FV. The signal kaon track 
reconstructed in the DC must satisfy the following requests: 



40 cm < < 150 cm. 



z\ < 150 cm 



(4.1) 



-RpCA < 10 cm, 

I^pcaI < 20 cm, 

70 < pK < 130 MeV. 



(4.2) 
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where -RpcA = y ^pcA + ^PCA PCA indicates the point of closest approach of the 
kaon track to the interaction point. The low level of background allows to keep these cuts 
very loose, even accounting for the effect of multiple scattering and energy loss of the kaon 
in the DC wall and into the beam pipe. We ask the reconstruction of a vertex formed 
by the kaon track and one of its charged decays, lying inside the FV. The only source 
of background is given by events in which one kaon track is split in two pieces by the 
reconstruction algorithm so as to mimic a decay vertex. This happens for kaon with low pt 
which describe cirles in the DC. In order to reject these events we used the momentum of 
the charged secondary particle {p*) evaluated in the kaon rest frame, using the kaon mass 
hypothesis. With the cut p* > 100 MeV we lose about 5.4% of signal and we reduce this 
background to the 0.46% level. Since charged kaons have an average velocity /3 ~ 0.2 they 
lose a non negligible fraction of their energy traversing the beam pipe and DC walls and 
in the DC gas. Therefore, given the decay length, we have to correct for the corresponding 
change in /?, of the order of 25%, to evaluate the proper time. Once the decay vertex has 
been identified, the kaon track is extrapolated backwards to the interaction point in 5 mm 
steps, Ali, taking into account the average ionization energy loss dE/dx to evaluate its 
velocity j3i in each step. The kaon proper decay time, t* is finally evaluated as 



i i 

In the range 12 < t* < 40 ns we collected 1.7 x 10^ events. The average resolution in t* , 
evaluated from MC simulation, is of the order of 1 ns, see fig. |l| and its effects are taken 
into account in the fit procedure by means of a resolution smearing matrix, as described 
in section |^, which makes the resolution folding. 

4.2 Kaon decay time method 

The second method relies on the measurement of the kaon decay time using EMC infor- 
mation only for the signal side. We consider events with a -k^ in the final state: 

^ X + 7r° ^ X + 77 (4.4) 

In this case the reconstruction efficiency of the kaon decay vertex is evaluated using a data 
control sample given by DC information only. We obtain the kaon decay time using the 
photon arrival time to the EMC. From the measured momenta of the meson and of the 
tagging kaon we can evaluate the momentum of the signal kaon at the IP and build the 
expected path of the signal kaon: it is obtained by geometrical step by step extrapolation, 
corrected for local magnetic field and energy loss effects. Then we look for clusters not 
associated to tracks in the DC, "neutral" clusters. Among these we select as photon 
candidates the two (or three, if present) most energetic clusters with: 

24° < 6'ci < 156° 

To < tci < To + 70 ns (4.5) 
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where Od is the polar angle, the angle between the position of the cluster with respect 
the IP and the z axis, td is the time of the cluster, and Tq is the decay time evaluated 
considering the JC^2 decay chain on the tagging side. These cuts are used to reject the 
machine background. The request of three neutral clusters allows us to select also the 
events with the track of the charged decay particle not associated to its calorimeter cluster. 
Then we move along the expected path of the signal kaon in 5 mm steps and at each step 
we look for the vr" — > 77 decay vertex by minimizing a function: 



In Eq. (4^) M^^ is the invariant mass of the two photon candidates, M^o is the vr^ 



mass and aM — ^''^ MeV is the resolution on the vr'^ mass; Ati2 is defined as: 

Ati2 = (^71 - n/c) - - r2/c) (4.7) 

where t^. is the time of the i-th neutral cluster, rj is the distance of the neutral cluster 
from the candidate kaon decay vertex; uj^j is obtained by propagating the errors on the 



equation 4.7. tx is the time of flight of the charged kaon measured along the kaon path, t^i 
is the time of flight of the charged kaon given by the weighted mean of the arrival time of 
the two photons; at is the uncertainty on the kaon decay time obtained from the selected 
neutral clusters, of the order of few hundreds ps. In case three clusters have been selected 
we choose the pair with best x^- The position along the signal kaon path that gives the 
minimum value of the defines the ^ X + vr'^ ^ X + 77 decay point and is accepted 
if lying inside the FV. Moreover we require: 

< 30, 

80 MeV < M^^ < 200 MeV, (4.8) 

Atl2 



< 5. 



The charged kaon proper time is obtained from: 

t* = ^ (4.9) 

(7) 

Where (7) is the average between the charged kaon Lorentz factor 7 at the IP and at the 
decay vertex. This approximation is appropriate because the variation of 7 along the decay 
path is of the order of 1%. The average t* resolution of the events selected is better than 1 
ns, see fig. |^. The secondary peaks at ibnx2.7 ns are due to events with an incorrect bunch 
crossing determination. In the range between 10 and 50 ns we collected 8 x 10^ events. 



5. Efficiency evaluation 

In order to be as little dependent as possible from the MonteCarlo simulation and from 
possible differences between data and MonteCarlo, we measured directly on data the recon- 
struction efficiency of the kaon decay vertex, as a function of the charged kaon proper time. 
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Figure 1: t* resolution from kaon path length measurement (left) and from kaon decay time measurement 
(right). The side lobes at ±nx2.7 ns in the right plot are events with incorrect bunch crossing time 
determination. This effect is absent for the length measurement. 



for both methods. A data control sample for the efficiency of the first method, using DC 
information for signal selection, is obtained from the signal selection of the second method, 
based on EMC information only, and viceversa. What is relevant for the fit is the efficiency 
behavior as a function of t* rather than its absolute value. In both cases the efficiency has 
been obtained in bins of t* using the definitions (O) and (O) respectively. We applied the 
same method on MonteCarlo events and we compared the efficiencies obtained with this 
method (MC datalike efficiencies) with the MC true efficiencies, determined as functions 
of the true proper time including all the decay modes used for each method [K^ — > all 
for the first method and — > Xtt^ for the second method). Fig. ^ shows a comparison 
between the MC true and the MC datalike efficiencies. For the first method it is important 
to stress that even if the control sample is obtained reconstructing kaon decays with a 
in the final state, datalike efficiency reproduces the true one obtained reconstructing all 
the kaon decay modes. Fig. |3| shows the ratio of the efficiency measured on data over the 
MC datalike efficiency, for both methods. For the kaon decay length method, the efficiency 
measured on data in the range [5, 20] ns is different in shape from the MC efficiency. The 
highly ionizing charged kaon ffi'es multiple hits in the small cells of the DC, i.e. in the inner 
12 layers of detector. This effect has been introduced in the MC simulation afterwards. 
For this reason we used the efficiency evaluated directly data with the small, O{10~^), 
correction given by the ratio of the MC true efficiencies over MC datalike efficiencies. 



6. Fit to the proper time distribution 

The fit to the proper time distribution is done comparing the observed with the expected 
distribution and minimizing a function. The entries in each bin of the expected his- 
togram are given by the integral of the exponential decay function (which depends only 
on the kaon lifetime) corrected for the efficiency curve; a smearing matrix accounts for the 
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Figure 2: Left: charged vertex reconstruction efficiency as function of the charged kaon proper 
time. Right: 7r° decay vertex efficiency as a function of the charged kaon proper time. Dots (open 
circles) represent MC datalike (MC true) efficiencies. 
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Figure 3: Left: ratio between the charged vertex reconstruction efficiency as function of the 
charged kaon proper time evaluated on data and the MC datalike one. Right: ratio between the 
tt" decay vertex efficiency as a function of the charged kaon proper time evaluated on data and the 
MC datalike one. 



effect of the resolution. We fit the proper time distribution in the time interval where we 
have good agreement between MC true and MC datalike efficiencies. The efficiency profiles 
obtained on data as in section |^ are corrected for the tiny residual difference between the 
MC true and MC datalike shapes, and for the small effect of the tag-signal hemispheres 
correlation; the overall corrections is of O(10~^) for both methods. The resolution function 
is evaluated in slices of the order of few nanoseconds on a MC data sample of about 175 
pb~^ of (/) — > K~^K~ events. We then define the N x N smearing matrix, whose element 
Sji gives the fraction of events generated in the i-th bin but reconstructed in the j — th bin. 
We chose a bin size of 1 ns in order to reduce the statistical fluctuations and the relative 
importance of the smearing corrections. We then minimize 

7=1 ^ i ' 



7 



where is the number of entries in the j — th bin of the proper decay time distribution. 
The content of the j — th bin of the expected histogram, (NJ^), is given by: 

N 

1=1 

where /^(t) is the integral of the exponential over bin i, is the measured reconstruction 
efficiency and €'^°^^ is the efficiency correction described above. The statistical fluctuation 
of the expected histogram in the bin j — th, Uj*, is given by the sum of the statistical 
fluctuation of the efficiency, the statistical fluctuation of its correction and the poissonian 
fluctuation of the expected histogram. Nf^i is the number of bins used in the fit. The 
number of bins of the expected histogram is larger than Nfn since we have allowed the 
migration from/to the bins used in the fit for a slightly larger range. 



6.1 Kaon decay length results 



We build the expected histogram in the region between 12 and 40 ns and we fit it in the 
region between 15 and 35 ns. For we obtained: 



r+ = (12.338 ± 0.042) ns 
While for the K~ we obtained: 

r" = (12.395 ± 0.045) ns 



xVndf = 18.2/18 P(x^) = 44.7% . (6.3) 



xVndf = 26.7/18 P(x^) = 10.8% . (6.4) 



The left panel of fig. ^ and |^ shows the data distribution of t* together with the fit results. 
The right panel shows the distribution of the residuals defined as the difference between the 
data distribution and the fit. For both charges the distribution of residuals is satisfactory 
fitted by a constant compatible with zero within errors. The measurements obtained for 
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the two charges are in agreement with each other. Their weighted mean is: 

T = (12.364 ±0.031) ns. 



(6.5) 
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Figure 5: Decay length method. Left: fit to the K~ proper time distribution: the black dots are 
data distribution, in grey the expected histogram fitting function. Right: Residuals in the region 
between 15 and 35 ns. 



6.2 Kaon decay time results 



For the second method we build the expected histogram in the region between 10 and 50 
ns. We made the fit in the region between 13 and 42 ns. The result from the fit to the 
proper decay time distribution, see fig. ^, is: 

r+ = (12.315 ± 0.042) ns xV^df = 25.5/27 P(x^) = 49%. (6.6) 

For the K~ lifetime, see fig. 0, we obtained: 

T" = (12.360 ± 0.043) ns xV^df = 22.9/27 P(x^) = 69%. (6.7) 



The measurements obtained for the two charges are in agreement with each other. Their 
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Figure 6: Decay time method. Left: fit to the proper time distribution: the black dots are 
data distribution, in grey the expected histogram fitting function. Right: Residuals in the region 
between 13 and 42 ns. 



weighted mean is: 



T = (12.337 ±0.030) ns. 



(6.8) 



As shown in the fig. ^, |5|, ^ and |^ the agreement between data distribution and expected 
histo distribution is very good also outside the fit range for both the charges and both the 
methods. 
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between 13 and 42 ns. 



7. Systematic uncertainties 

All the measurement strategy has been designed in order to have redundancy to keep 
systematic effects under control, and also to use as much as possible control samples selected 
on data, rather than relying on MonteCarlo simulations, for all the inputs to the fit. The 
determination of efficiencies has been already discussed; the efficiency correction effect has 
been checked by performing the fit with and without it. For the resolution functions, 
which at first order we obtain from MonteCarlo, we have checked the data/MC agreement 
on the distribution of the difference between the t* from the decay length and the t* from 
the decay time, on doubly reconstructed events. This allows us to check also the amount 
of incorrect bunch crossing determinations since the length is clearly not affected by this 
feature. The core of the time resolution is also kept under control by comparing the time 
of the two neutral clusters. All checks have shown very small data/MC corrections and 
negligible variations to the fit results and quality. The stability of the fit procedure with 
respect to the fit range has been checked by varying significantly the fit window and by 
changing the bin size from 1 ns to 0.5 or 2 ns. The systematics due to the selection have 
been obtained by comparing the fit results with different cuts for both the signal and the 
control samples. The beam pipe and DC wall thickness are known to about 10%; this in 
turn affects the evaluation of the kaon energy loss in the materials and consequently the 
kaon /?. We have repeated both measurements considering different thicknesses (by varying 
them within their errors) in order to assess the systematic effect on the lifetime. The energy 
loss inside the DC is more than one order of magnitude smaller than the one in the walls 
and the uncertainty on the material crossed inside the DC gives a negligible contribution 
to the systematic error on the lifetime measurements. The first method is very sensitive 
to systematic shift on the (3 of the kaon (see eq. [4. 3D and we observe a sizeable effect on 
the lifetime. This is not true for the second method (see eq. |4.9| ): for a typical charged 
kaon in our detector 5'y/j ~ 0.06 x 5P/P and td depends on the kaon energy loss only at 



second order, via the third term of the function in eq. 4^. For the time method the 
overall systematic effect due to the uncertainty on the kaon energy loss has been actually 
found negligible. The estimates of the systematic effects are summarized in table |l|. The 
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systematics of the two measurements are almost completely uncorrelated. The result from 



systematic uncertainty 


Length (ps) 


Time (ps) 


Fit range 


± 10 


± 10 


binning 


± 10 


± 10 


efficiency correction 


± 15 


± 10 


beam pipe thickness 


± 10 


negligible 


DC wall thickness 


± 15 


negligible 


Signal sclecticni cuts 


± 15 


uogligiV)lc 


Control sample cuts 


negligible 


±10 


Resolution effects 


negligible 


negligible 


total systematic uncertainty 


± 31 


± 20 



Table 1: Sources of systematic uncertainties, 
the length measurement is: 

r = (12.364 ± 0.031stat ± 0.031syst) ns 
and from the time measurement is: 

r = (12.337 ± 0.030stat ± 0.020syst) ns 

8. Correlation and average 



(7.1) 
(7.2) 



In order to average the two methods we calculate the statistical correlation between the 
two results. The correlation arises because the data samples used for the two methods 

have ~ 30% of events in common. The normalized correlation is 30.7% in agreement with 
a direct estimate obtained dividing the data in subsamples. Assuming that the systematic 
uncertainties are uncorrelated we then obtain the average: 



T = (12.347 ± 0.030) ns 



(8.1) 



The measurement obtained agrees, within the errors, with the result given by Ott and 
Pritchard, [3]: 

r = (12.380 ± 0.016) ns (8.2) 

and with the PDG fit, [4]: 

T = (12.384 ± 0.024) ns (8.3) 

9. CPT test 

The comparison of and lifetimes is a test of CPT invariance which requires the 
equality of the decay lifetimes for particle and antiparticle. The average of the two methods 
is: 



r+ = (12.325 ±0.038) ns 
r" = (12.374 ± 0.040) ns 



(9.1) 
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From these measurements, and taking into account that most of the systematic effects 
cancel out in the ratio, we obtain: 

^ = 1.004 ± 0.004 (9.2) 

Our result agrees with CPT invariance at the 0.4% level. Ref [10] had already verified 
agreement at the 0.08% level. 
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